Genuine multipartite nonlocality of entangled thermal states 
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We assess quantum non-locality of multiparty entangled thermal states by studying, quantitatively, both tri- 
partite and quadripartite states belonging to the Greenberger-Horne-Zeilinger (GHZ), W and linear cluster-state 
classes and showing violation of relevant Bell-like inequalities. We discuss the conditions for maximizing 
the degree of violation against the local thermal character of the states and the inefficiency of the detection 
apparatuses. We demonstrate that such classes of multipartite entangled states can be made to last quite signifi- 
cantly, notwithstanding adverse operating conditions. This opens up the possibility for coherent exploitation of 
multipartite quantum channels made out of entangled thermal states. Our study is accompanied by a detailed 
description of possible generation schemes for the states analyzed. 
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I. INTRODUCTION 

Where lays the boundary between classical and quantum 
worlds? The daunting nature of this question is not preventing 
the pursuit of very interesting studies digging deeply into the 
origin of non-classicality of a physical system. On the con- 
trary, the difficulties inherent in such a fundamental investiga- 
tion are sharp stimuli to the research of tests and physical sys- 
tems able to challenge the common belief that quantumness 
occurs only under quite special, and yet hard to achieve, con- 
ditions. Very recently, quite significant endeavors have been 
produced in the study of quantumness at the "large scale" by 
proposing ways to infer non-classicality in multi -photon states 
and massive mechanical oscillators (QJ-Q] . 

In this article we contribute to such a quest by address- 
ing the noteworthy case of multipartite quantum correlations 
shared by systems that, when individually taken, are fully 
classical and revealed by instruments far from offering any 
single-quantum resolution. We address the case of multipar- 
tite entanglement shared by bosonic systems that are locally 
prepared in chaotic thermal states, which are commonly in- 
tended as well-defined classical entities, and dub them en- 
tangled thermal states (ETS) [8]. We demonstrate that such 
quantum correlations can be easily made strong enough to be 
revealed, through the violation of suitable Bell-like inequal- 
ities, against any initial local temperature and regardless of 
the coarse-graininess of the detectors used in order to imple- 
ment the non-locality test. We first study ETS versions of tri- 
partite GHZ [19] and W states [ 10], which are prominent and 
non-equivalent classes of three-particle entangled states. We 
discuss two schemes for generating GHZ-like ETS and draw 
a comparison between the slightly different states therefore 
obtained. We show that an inequality for genuine multipartite 
non-locality can be violated up to the maximum value allowed 
for a given representative of the class of GHZ-like ETS, when 
effective local rotations and highly inefficient and noisy ho- 
modyne measurements are employed. Similar action can be 
taken for the ETS version of W states llfJI . where the viola- 
tion is not to the maximum degree. We then extend our anal- 
ysis to larger entangled states, proving independence of a few 



of our results from the number of bosonic systems involved in 
the multipartite states we scrutinize. Such larger ETS states 
include the quadripartite version of GHZ-like ETS and a four- 
mode linear cluster-like ETS. Our results go along and are 
consistent with a recently-started line of investigation aiming 
at showing that quantumness can be enforced and revealed in 
situations that are at the verge of the classical world lf5ll8i [Tlll . 

The remainder of this paper is organized as follows. In 
Sec. [II] we assess genuine tripartite non-locality in GHZ and 
W-like ETS. The main tool of our investigation is an inequal- 
ity derived by Svetlichny IU2I1 . We sketch two schemes for 
the generation of GHZ-like ETS states and provide an opera- 
tive protocol for testing tripartite non-locality by means of ef- 
fective rotations and arbitrarily efficient homodyne measure- 
ments. In Sec. |III]we tackle the non-local properties of four- 
party states such as quadripartite GHZ-like ETS and the in- 
teresting class of linear cluster-like ETS. Non-locality can be 
revealed regardless of the local temperature and can be made 
robust to the inefficiency and fuzziness of the detection de- 
vices. Finally, Sec. HVIsumm arizes our findings. 



II. TRIPARTITE CASE: OF GHZ AND W-LIKE ETS 

A. Svetlichny versus Mermin inequality for spin- 1/2 particles 

The quantification of the degree of genuine multipartite 
entanglement in the aforementioned class of tripartite and 
quadripartite ETS is quite a daunting task. For pure states of 
three entangled qubits, it is possible to quantify the degree of 
multipartite entanglement jl3[]- However, this is not the case 
for the class of ETS here at hand. In fact, as it will be seen 
in Sec. Ill Bl we will be dealing with highly mixed states of 
systems spanning infinite dimensional Hilbert spaces. More- 
over, by construction, such states do not belong to the class 
of so-called Gaussian states, i.e. states whose Wigner func- 
tion is a Gaussian in phase-space lfl4ll . While correlations in 
Gaussian states are well and easily characterized, we face the 
lack of necessary and sufficient criteria for the quantification 
of entanglement in non-Gaussian states. In fact, the available 
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FIG. 1: (a) Conditional scheme for generating GHZ-like ETS. Three displaced thermal states (labelled j = 1, 2, 3) interact with a two-level 
ancilla prepared in state \+) A = (|0) A + |l) A )/\/2. The interaction is ruled by the coupling Hamiltonian Haj = hQ.\l) A(l|ajaj for a 
time equal to n/Q,. The ancilla is finally projected onto |±} A with a(+|— )a = 0. (b) Beam splitter-based scheme for the generation of 
tripartite GHZ-like ETS. A displaced thermal state as introduced in the body of the paper undergoes a transformation that puts it into the form 
p' h '(V,d) = A/+ / d 2 aP l a(V, d) \s+) (s+\ with |s+) defined as in Eq. $4$ and A/+ a normalization constant. This is then superimposed to 
two vacuum modes at beam splitters of transmittivity t\ — a/2/3 and t2 — l/\/2 1 15]. (c) Scheme for Svetlichny test performed over a 
GHZ/W-like ETS. Each mode of a state prepared off-line is first appropriately rotated and then projected onto in-phase quadrature eigenstates 
by means of arbitrarily efficient homodyne detectors. The outcomes of the measurements are appropriately dichotomized (see the other inset). 



entanglement measures in continuous variable are based (to 
the best of our knowledge) on the use of the negativity of 
partial transposition criterion formulated in terms of covari- 
ance matrices, which carry exact information on the state of 
a system only in the Gaussian scenario. Therefore, even the 
case of quantifying bipartite entanglement represents quite a 
challenging problem. A partial assessment of the entangle- 
ment content of a given two-mode state belonging to such 
class passes through so-called entangling power methods: one 
could establish a lower bound to the entanglement within a 
two-mode non-Gaussian mixed state by determining the en- 
tanglement that can be transferred to two qubits by means of 
bi-local unitaries involving a mode and a qubit. Such a strat- 
egy has been used in Ref. [16], where the entanglement in a 
mixed two-mode non-Gaussian state has been examined un- 
der the assumption that each mode-qubit subsystem interact 
via a local Jaynes-Cummings model. Although the technique 
is fundamentally interesting in light of the intricate relation- 
ship between systems of different dimensions, it is of not easy 
generalization to the multipartite scenario and we can at best 
provide only a lower bound to non-classical correlations, an 
exact quantification still being quite elusive. The lack of effi- 
cient ways to quantify entanglement motivates an assessment 
of non-classicality of ETS in terms of multipartite non-locality 
very appealing and, basically, the only rigorous way to ascer- 
tain whether or not a given non-Gaussian state of many modes 
is righteously entangled or not. This is the perspective which 
we will take in our study. 

Although one may be tempted to identify one with the 
other, quantum non-locality and quantum entanglement in 
multipartite settings are concepts which should be approached 
carefully. In fact, it is straightforward to realize that the vio- 
lation of an n-particle Bell's inequality by an n-particle en- 
tangled state is not sufficient to guarantee that genuine multi- 
partite entanglement is shared by the system's elements. The 



non-local nature witnessed by the violation of a Bell's inequal- 
ity might well be entailed simply by an entangled state involv- 
ing only m < n particles. In 1987, Svetlichny addressed this 
point by deriving an inequality, for the tripartite case, that is 
able to discriminate between two- and three-particle entan- 
glement lfl2ll . Such inequality is satisfied by restricted local 
realistic models allowing for a degree of two-particles non- 
locality. 

We now give a brief account of the inequality formulated 
by Svetlichny. Let us consider Mermin's version of Bell's in- 
equality for three spin- 1/2 particles [ 17]. In one of its possible 
formulations, Mermin function A4 is built from the correlation 
function C(abc) as 

M = C(abd) + C{ab lC ) + C{aibc) - C(aihci), (1) 

where the pairs (fe, k\) with k — a,b,c are the two dichotomic 
outcomes of an observable measured at the detection stage for 
one of the particles. Here, C(abc) is the statistical correla- 
tion function for measurements having outcomes a, b and c 
respectively. For local realistic theories (LRT), it is < 2. 
By exchanging k with k\ in Eq. ((TJ, one gets a new Mermin 
function A4 i that, under LRT assumptions, obviously satisfies 
(in modulus) the same bound as A4. Therefore, by taking the 
Svetlichny function S (a, b, c) = \A4 + we get 

5(a,b,c) = \C(abcx) +C(ab lC ) +C{a 1 bc) +C(abc) 
— C(a 1 6 1 c) — C(aibci) — C(abiCi) — C(a 1 & 1 c 1 )| < 4, 

(2) 

where k = (fc, fci ). Quantum mechanics, on the other hand, 
predicts the existence of genuinely tripartite quantum cor- 
related states violating such a bound. In particular, when 
the correlations C(abc) are evaluated over a tripartite GHZ 
state H, the value 4\/2 is obtained for the Svetlichny func- 
tion, which is the maximum value achievable for any tripar- 
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tite state. Svetlichny inequality (SI) is maximally violated by 
projecting each particle j in a GHZ state onto the eigenstates 
|±) ■ of the observable cos •do x + sin day — 9,fap), where 
a x,y,z are the three Pauli matrices. A straightforward calcula- 
tion gives the correlation function 

C(0,fap) = cos(6» + + p). (3) 

For 6 = 6»i + tt/2 = 37r/4, 4> = -/x x = tt/2 and p = fa = 0, 
we have \M\ = \Mi\ = 2y/2 and \S\ = 4^2, which shows 
violation of both Svetlichny and Mermin inequality. As dis- 
cussed by Cereceda fHHl . SI is a righteous Bell inequality for 
the tripartite case and emerges as a valuable tool for the unam- 
biguous ascertain of the existence of genuine tripartite entan- 
glement and tripartite non-locality for any three-particle state. 
To this task, the use of a standard Mermin's inequality IU7I1 is 
not sufficient: on the contrary, quantum correlations that vio- 
late SI are strong enough to maximally violate Mermin's in- 
equality as well. The inequality by Svetlichny has been inde- 
pendently extended to the n-partite scenario in Refs. lfl9ll20ll . 
Very recently, Lavoie et al. have experimentally demonstrated 
the violation of SI by a tripartite GHZ state l^] encoded in the 
polarization degrees of freedom of three photons in a linear 
optics setup 12 ill . 



B. Generation of GHZ/W-like ETS and violation of Svetlichny 
inequality 

We now address tripartite non-locality of ETS. In our pro- 
posal, three local parties are each provided with one mode 
j = 1, 2, 3 of a tripartite ETS that has been off-line prepared 
using single-mode displaced thermal states ||8t ll ill . These 
are defined as pf(V,d) = J d 2 aP^ h (V,d)\a} j (a\, where 

P^ h (V,d) = 2[n(V - l)]-^- 1 ^^ is a Gaussian func- 
tion with center d (with respect to the origin of the phase 
space) and variance proportional to V = 2n + 1, where n 
is the mean thermal occupation number of the mode. Here, 
\a)j is a coherent state of mode j, which has associated 

bosonic operators dj and a-. Displaced thermal states are 
the building blocks for the construction of a tripartite GHZ- 
like ETS state, as described in Figs. Q] (a) and (b). The 
scheme in panel (a) is probabilistic and based on condi- 
tional displacements of each of the three modes upon inter- 
action with a two-level ancilla A having logical states |0) A 
and \1) A . This is realized by enforcing the mode-ancilla 
coupling Haj — and upon preparation of 

A in \+) A = (\0) A + \l ) A ) '/y/2. Nonlinear media with 
free-traveling optical fields 11221 |23J or dispersive interactions 
within optical/microwave cavities J24JL. ma Y be used to im- 
plement the required interactions flU, [la] . The state of A is 
eventually projected onto the state basis {|+}^ , |— ) A } (with 
a{~ \+)a = 0), as shown in Fig.[T](a). The scheme in Fig.Q] 
(b), on the other hand, relies on interference at appropriately- 
arranged beam splitters, in a way completely analogous to the 
proposal put forward in Ref. 1151 . Having off-line prepared 
state p th ' (V,d) = 7V+ / d 2 aP^(V,d) \s+) (s+\ of mode 1 



with 

|*±)oc(|a)±|-a», (4) 

a three-mode GHZ-like ETS is achieved by superimposing it 
to two additional modes in their vacuum state. State p th (V, d) 
can be prepared similarly to what is shown in panel (a), i.e. 
by letting a displaced thermal state of mode 1 interact, ac- 
cording to Hai and for a time 7r/fi, with a two-level ancilla 
prepared in \+) A - The latter is eventually projected onto the 
basis {\+) A , \-) A } to leave mode 1 in p th '(V,d) [8]. The 
scheme in Fig.[T](b) creates the state 

P ( gLi=^ J d 2 aP l *(V,d)\GHZ(a,a,&)) 123 (GHZ(&,a,a)\ 

(5) 

with a = a/v3, while the conditional approach in Fig.Q](a) 

gives 

P%,2 =^2 J d 2 ad 2 t3d 2 (P*(V,d)P$(V,d)P*(V,d) 
x |GHZ(a,/3,C)) 123 (GHZ(a,/3,C)| 

(6) 

with Wi.2 being normalization factors and a,/3,^€ C. In 
these equations we have used the (unnormalized) GHZ-like 
entangled coherent state 

|GHZ(a, P, 0) 123 = (1 + 2>) =1 e^ )\a,0, C) 123 (7) 

where 1 is the identity operator. Although conceptually analo- 
gous, the cases encompassed by Eq. © and (O are technically 

(3) 

different and the class of states represented by p K gf ^ z 2 turns out 
to be analytically easier to treat. Therefore, in order to provide 
a clear interpretation of our results, in what follows we study 
the case of Eq. (O in detail, leaving the assessment of Eq. §5§ 
to the numerical results shown in Fig. [3] It is worth noting that 
both the tripartite GHZ- and W-like entangled coherent states 
have been shown to violate Mermin's inequality llal . 

Our task here is to give a clear-cut account of the main re- 
sults of our investigation, providing at the same time an intu- 
ition of the physical features behind the observation of non- 
classical multipartite correlations in the classes of multiparty 
ETS-based states addressed here. We aim at showing violation 
of SI under unfavorable conditions such as high-temperature 
local bosons and "fuzzy" measurements. We address the case 
of a set of local rotations performed over each of the modes 
entering a GHZ-like ETS state. Explicitly, we consider the set 
of angles (8, 7) and the associated transformation 

^ / sin(0/2) e^cos(9/2)\ 

Rl{9 ^>-{e-^co S (d/2) -sin(0/2) ) (8) 

to be applied to the vector of coherent states (\a) 1 \~a)- i ) T 
for mode 1. Modes 2 and 3 experience similar transforma- 
tions, each determined by the pairs of angles (fa 5) and (p, v) 
respectively. An implementation of such a set of rotations 
using non-linear interactions and phase-space displacements 
is discussed in Refs. iTTTL l25ll Following the arguments by 
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Svetlichny [12] illustrated above, we restrict the set of local 
rotations over the modes at hand to6> = = /j = 7r/2 and, un- 
less stated otherwise, indicate the rotation operators as Rj (a) 
(a = 7, 8, v). The projections needed in order to evaluate joint 
probabilities and correlation functions as described in Eq. (O 
are instead implemented by dichotomizing the outcomes of 
homodyne measurements performed over the three modes in 

(3) 

Pghz 2- W e associate a logical +1 (—1) to a positive (negative) 
homodyne signal, as illustrated in Fig.Q](c). For the tripartite 
case, we calculate the correlation function as 



C( 7 ,<5, v) = R 



+++ 



P___+II(.P + __ (9) 



where II represents a permutation of the pedices appearing 
in the conditional probabilities Pm p (with k,l,p = ±). The 
latter are easily calculated as 



Pklp 



dx dy dz {x,y,z\pf^ z2 \x,y,z), (10) 
iz k Jn, Jn TJ 



where we have introduced the rotated GHZ-like ETS 

P%,2 = ^(T^WflsM&All^^M and the 
positive (negative) domain of integration 1Z+ = [0,oo) 
(1Z- = (— oo, 0]). It is thus found that 

Cghzflil, v) = Nghxpil, 6, v) cos(7 + S + v)Erf 3 [^£] 

yv 
(11) 

with Erf[x] the error function of argument x. The appear- 
ance of the error function in the result above is noteworthy. 
It arises from the finite-range integration over the products of 
the Gaussian P* function entering the definition of a thermal 
state and the one resulting from the homodyne probing of the 
fields' quadratures. When d > y/V/2, Erf[V2d/W] -> 1 
with Afghz.2 —> 1> making the correlation function identical 
to the one obtained for the case of a tripartite spin- 1/2 GHZ 
state 112 111 . This is the key of our results. As confirmed by the 
analysis performed in the remainder of the paper, we can gen- 
eralize the functional form in Eq. (fTTb and claim that the cor- 
relation function for any ETS-based multipartite state will en- 
compass an angular part that, under ideal rotations, is identical 
to the one achieved for spin-1/2 particles. The modulus of the 
correlation function, though, also depends in a crucial way on 
the degree of distinguishability of the multipartite-state com- 
ponents, as encompassed by the error function. This accounts 
for the clear universal behaviour of the correlation functions 
associated with ETS-based states against the key parameter of 
the local thermal components. 

When inefficient homodyne detectors are considered, all 
having the same detection efficiency r), the correlation func- 
tion is easily found from C g h z ,2 (l, 8, v) with the replacements 
d — > rjd and V — > 1 + r] 2 (V — 1) in the error function appear- 
ing in Eq. ( fTTb . So, the effect of an inefficient detector is just 
to slow down the saturation of the ^-dependent term in the 
correlation function. This simply means that, for set values of 
V and for 77 < 1, a larger value of d is required to achieve 
the maximum allowed degree of violation of Svetlichny in- 
equality. More in detail, for ^ > 1 we need to displace 



the local thermal states by d v ~ dy/ 1 + (Vr/ 2 )^ 1 in order 
for the corresponding correlation function to reach the value 
of Cghz^il, 8, v) corresponding to 77 = 1 (i.e. for perfectly 
efficient detectors). The scaling is quite favourable at large 
values of V, which shows that highly thermal states that are 
initially considerably displaced with respect to the origin of 
phase space are extremely insensitive to the effects of ineffi- 
cient homodyne detection. One could even ignore the actual 
detection inefficiency and, at large V, displace the local states 
as if the homodyners were ideal, without affecting the perfor- 
mance of the scheme. From the above discussion it is straight- 
forward to understand that, under the present conditions and 
for the GHZ-like ETS state addressed above, we get 



S« h z < 2 (7, 8, u) oc S ghz ' 2 (-y, 8, iv)Erf 3 [- 



V2di] 



yi + ^-i) 

(12) 

where S ghz,2 (') 1 8, v) is the Svetlichny function for the ideal 
case of spin-1/2 particles and the proportionality sign is due 
to a factor that tends to 1 as d 2 > (r)~ 2 + V - l)/2. This 
shows that non-classicality as witnessed by the violation of 
Svetlichny inequality will be observed up to its maximum al- 
lowed degree even for large initial temperatures of the mul- 
tipartite state. Fig. [2] shows the Svetlichny function against 
V and d for 77 = 0.1. Maximum violations of the inequal- 
ity can be observed for any value of V by choosing a suf- 
ficiently large d. The parameter d is the "knob" to tune in 
order to optimize the non-classical properties of a given state 
at an assigned value of the thermal spread V. Intuitively this 
means that the entanglement in an ETS is a delicate trade- 
off between the distinguishability of its state components, as 
measured by their mutual distance d in phase-space, and the 
width V of each thermal distribution. When the Gaussian 
probability functions defining each thermal state are so large 
that they overlap significantly, the state components become 
quasi-indistinguishable and entanglement is correspondingly 
destroyed. It should thus be clear that, per assigned value of 
V, a way to counteract such entanglement washing-out effect 
is to make the state components sufficiently distinguishable in 
phase space, which implies the increase of d. 

An analogous behavior is found also for the alternative form 

(3) 

of GHZ-like ETS pV r In this case, though, due to the 
complications related to the calculations that are necessary to 
get the correlation functions, an analytic expression for the 
Svetlichny function turns out to be not possible. Neverthe- 
less, the inequality can be well assessed, quantitatively, by 
relying on a numerical approach that allows us to assess the 
behavior of S^ hz,1 (j, S, v) against d, at set values of V and 
77. While, qualitatively, the very same features highlighted for 
the case of 2 hold also in this case (monotonic growth 
of the Svetlichny function with d and saturation to the degree 
of maximum violation for d 3> W), it is interesting to de- 
termine if one of the state-generation strategies proposed here 
offers any advantage when evaluated against the violation of 
Svetlichny inequality. In Fig. |3]we show a quantitative com- 
parison between S ghz,1 (-f, 8, u) and S^ hz ' 2 (j, 8, v) for two 
values of V (arbitrary choices) and using the same detection 
inefficiency 77. We can clearly see that, despite the evident 
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having the form p th "{V,d) oc J d 2 aP^{V,d)\r+)(r+\ is 
achieved, where |r+) oc (\a)+i \— a)). We have checked that, 
for the GHZ-like ETS that is generated by using p th (V, d) in 
the protocol of Fig. [TJ (b), there is always a set of local ro- 
tations that, complemented by dichotomized homodyne mea- 
surements, allow for the violation of the tripartite Svetlichny 
inequality in a way fully analogous to what has been ad- 
dressed so far. 



FIG. 2: Violation of Svetlichny inequality by a tripartite GHZ-like 
ETS under effective rotations and inefficient homodyne detection (ef- 
ficiency 77 = 0.1). We show the Svetlichny function against displace- 
ment d and variance V of the local thermal distributions. The floor 
of the plot is given by the local realistic bound of 4. For d ^> vV, 
the upper bound of 4\/2 is achieved, exactly as it would be for a pure 
GHZ state of three spin-1/2 particles under sharp measurements. 



similarities in the form of the two Svetlichny functions, the 
generation protocol illustrated in Fig. Q] (b) turns out to be 
slightly more convenient: the local realistic bound of 4 [see 
Eq. (01] is surpassed for slightly smaller values of d, when the 
beam splitter-based generation scheme is used. Such a trend is 
shown regardless of the value of V > 1 and 77 chosen for our 
numerical test and we have strong numerical evidence that the 
distance between the Svetlichny functions associated with the 
two generation schemes opens up (at small values of d) as V 
increases. In Sec.[]Tj]we will see that a similar result holds for 
a quadripartite linear cluster state. It is worth stressing the ex- 
istence of another possibility to generate tripartite GHZ-like 
ETS by means of the scheme in Fig. [TJ (b). Instead of re- 
lying on a resource having the form of p th (V, d), one could 
well use state p th (V, d) = Up lh (V, d)!/' 1 , where the unitary 
evolution U = e~ Z7T< - a a ) 2 / 2 is obtained by means of a self- 
Kerr non-linear medium. In this case, a single-mode resource 




FIG. 3: Comparison between the Svetlichny function associated with 

(3) (3) 

state Pgji z 1 (full lines) and p g ^ z 2 (dashed ones) for V = 5, 10 and 
rj = 0.3 (assumed to be the same for every detector). 
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FIG. 4: Svetlichny function for a tripartite W-like ETS at various 
values of V. The local realistic bound is 4 (shaded region of the plot). 
The maximum achieved violation of Svetlichny inequality, reached 
at large values of d, is quantitatively the same as in the spin-1/2 case. 



The qualitative features discussed above do not depend on 
the form of the multipartite state being considered. In fact, 
by following arguments similar to those valid for a GHZ-like 
ETS state, one can easily verify that the correlation function 
corresponding to an ETS version of the W state, which can be 
generated as discussed in lfl5ll . reads (for 77 = 1) 

, V2d 

C w (8, 0, p)oc [cos#cos0cos/j+2cos(# + </>+/u)]Erf [ — ]. 

VV 
(13) 

The angular part is identical to what has been ob- 
tained by Cereceda [18] for a tripartite spin-1/2 W state 
(l/v3) Z~2j=i &x,j |000) 123 - This class of states is known 
to violate Svetlichny inequality by ~ 0.355 lUsll . which 
is achieved by projecting each party onto the eigenstates 
of cos "!9o- 2 + sin-^iTa; (•& — 6,<j),n). In our formal- 
ism, this is equivalent to applying the local rotations 
Rj (2 arctan[(cos d — 1)/ sin d] , 0) to each mode and perform 
dichotomic homodyne measurements as described above. The 
corresponding Svetlichny function Sw is shown in Fig.|4]for 
several values of V and 77 = 1. Evidently, a large enough 
ratio d/VV makes a W-like ETS state violate the tripartite 
Svetlichny inequality up to the maximum degree allowed to 
the spin-1/2 counterpart of such states. The inclusion of de- 
tection inefficiency and the corresponding results are perfectly 
analogous to what was discussed in this Section regarding the 
GHZ-like ETS case. 
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in. QUADRIPARTITE CASE: GHZ AND CLUSTER-LIKE 

ETS 

A. Quadripartite GHZ-Iike ETS 

We now extend our study to the case of quadripartite ETS. 
As a first, simple step, we address the extension of our previ- 
ous study to non-locality of a quadripartite GHZ-like ETS. As 
done in the tripartite case, we focus our attention to the general 
state-engineering scheme depicted in Fig.[T](a), which can be 
straightforwardly generalized to the four-partite scenario by 
considering a fourth local displaced thermal state and, corre- 
spondingly, an additional ancilla-mode interaction. For equal 
local temperatures and displacements, the resulting quadripar- 
tite GHZ-like ETS takes the form 

p£l=^3 1 d 2 ad 2 pd 2 1 d 2 6P^(v,d)p^(v,d)p^(v,d) 

x Pf{V, d) |GHZ(a, fi, 7, 8)) 1234 (GHZ(a, j 9, 7, 5) \ 

(14) 

with A3 a normalisation factor, a, j3, 7, 6 € C and 




s\ 



FIG. 5: Svetlichny function for a four-partite GHZ-like ETS. As be- 
fore, we have used rj — 0.1. Local realisitic theories enforce a bound 
of 8 while the maximum achievable quantum mechanically is 8\/2. 

|GHZ(a,/?,<5, 7 )) = |a,/3, 7 ,5) + \-a, -7, S). In 
Ref. IU9I1 . Seevinck and Svetlichny have generalized the orig- 
inal argument by Svetlichny and provided the explicit form 
for the inequality to assess in the quadripartite scenario. This 
reads 

5(a, b, c, d)=\C(abcd) - C{a 1 bcd) - C(ab lC d) ~ C{abcid) 
+C(aibicdi) + C(aibcidi) + C{ab\C\di) + C(aibiC\di) 
—C(abcdi) — C(aib\cd) — C{a\bcid) — C{aibcd\) 
—C{abiC\d) — C(ab\cdi) — C{abc\di) + C{aib\Cid)\ < 8. 

(15) 

Quantum mechanically, there are genuinely multipartite en- 
tangled states that violate the local realistic bound of 8. In 
particular, a four-party GHZ state violates this inequality max- 
imally, achieving \S(a, b, c, d)| = 8^2 for the same set of 
local operations considered in Sec.HU When tested using the 
state in Eq. (14\ , the four-party Svetlichny function shows a 
behavior similar to what is shown in Fig. [2] for any speci- 



fied value of V and of the detection efficiency ?y, a sufficiently 
large displacement ensures the violation of SI up to the maxi- 
mum allowed value 8"\/2- Fig. [5] shows a significant example 
of such a behavior. 



B. Quadripartite linear cluster-like ETS 

Unlike the tripartite case, where only two inequivalent 
classes of multipartite entangled states exist [26], a plethora 
of choices are available in the quadripartite scenario. Among 
the states that can be considered, the class of cluster states is 
dressed of considerable interest in virtue of the role they play 
in the one-way paradigm for quantum computing [27]. It is 
thus interesting and relevant to assess the case embodied by 
such states and to build up their ETS counterpart. In what fol- 
lows, we focus our attention to the analogue of the spin- 1/2 
cluster state 

|C 1/2 ) = (|0000) + |1100) + 10011) - |llll»/2 (16) 

with {|0) , |1)} the logical states of each particle and the sub- 
script 1/2 reminding us that this form holds for a register of 
spins. Schemes for the generation of coherent state-encoded 
linear cluster states having a structure similar to the above 
have been put forward [28] and their usefulness as quan- 
tum channels for teleportation has been assessed by one of 
us l29ll . Here we describe suitable generation protocols for 
linear cluster-like ETS and then use them to test genuine 
quadripartite non-locality. The first scheme we propose is 
a straightforward modification of both the protocol put for- 
ward in Ref. ll28"ll and the general procedure utilized in or- 
der to generate GHZ-like ETS [cfr. Fig.[T](a)]. It is based 
on conditional displacement operations actuated by two ancil- 
lary two-level systems resulting from the same sort of interac- 
tion Hamiltonian introduced in Sec. [TTJ Differently from the 
scheme for GHZ-like ETS, here a mutual coupling between 
the ancillary systems has to be arranged according to H12 = 
^111)^45(111, where A and B are labels identifying the two 
ancillae, whose state after the network of couplings depicted 
in Fig.[6](a) is projected onto the {|+) , |— )}a,b bases. Upon 
detection of |++)^, the state of modes j — 1,2,3 and 4 
becomes 

Pauster= J d 2 ad^d^d 2 5P^(V,d)P^(V,d) 

x P*(V, d)Pf{V, d) \C(a, P, 7 , 5)} 1234 (C(a, 0, 7, S)\ 

(17) 

with a,j9, 7, 5 e C. By identifying a logical |0) (|1)) state 
of mode 1 with a coherent state of amplitude \a) 1 (|— a) 1 ) 
and using an analogous encoding for the remaining modes, 
state |C(a, j3, 7, 6)) turns out to be the coherent state-encoded 
linear cluster state 

|C(a, (3, 7, -5)) = i(|a, P, 7, S) + \a, (3, -7, -6) 
+ \-a,-j3,i,5) - \- a ,-p,-i,-S)). 
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FIG. 6: (a) Conditional scheme for the generation of a four-partite linear cluster state via local thermal states, conditional displacements and 
postselection. The interpretation of the symbols is the same as in the caption of Fig. Q] Differently from the GHZ-like ETS case, here we 
need two ancillary two-level systems and their mutual interaction. In panel (b) we show an alternative way of generating a linear cluster-ETS 
(see Ba an and Hoa in Ref. 1 34]). Two thermal states of annihilation operators &j (j = 1, 2) enter a non-linear cross-Kerr medium where the 
interaction H x (3) = hxa\<i\a\a2 takes place for a time n/x- The output modes are then superimposed to two ancillary modes prepared in 
their vacuum state at 50 : 50 beam spitters. Local phase shifters (not shown) complete the generation process. 



We will address an alternative state-engineering procedure 
in the next Subsection. The non-locality properties of spin- 
1/2 cluster states have been addressed in a seminal work by 
Scarani, Acin, Schenck and Aspelmeyer [30], which we dub 
from now on as SASA. Motivated by the fact that a linear 
cluster state does not maximally violate a four-partite Mermin 
inequality, SASA looked for an inequality that is maximally 
violated by quadripartite linear cluster states. By using the 
stabilizer operators, they have formulated a simple four-term 
correlator, here named SASA function. More in detail, for the 
form of the linear cluster state under scrutiny here, the spin- 
1/2 SASA function is obtained as the expectation value of the 
correlation operator 



0SASA = O- Z ,1 ® 



12 ® 0"x,3 



09 <7 y ,3 <&> Cj/,4 

(19) 

calculated over the cluster state. Classically, by associat- 
ing ±1 to any of the operators involved in Osasa, it is 
straightforward to check that the inequality |(Osasa)| < 2 
holds. Quantum mechanically, by construction, we get 
(C1/2 1 Osasa IC1/2) — 4. No other state can achieve a larger 
value, thus making the SASA inequality optimal for the class 
of states under scrutiny. Remarkably, such inequality has been 
experimentally tested using the polarization degree of free- 
dom of four photonic modes 113 ill . 



Our aim now is to prove that state p c i us ter defined in 
Eq. (T% violates SASA inequality up to its maximum value 
when the local thermal states are sufficiently displaced in 
phase space with respect to the variance V. As done before, 
the necessary correlators are implemented by dichotomising 
the outcomes of homodyne measurements performed over the 
four modes j = 1, 2, 3, 4 and using appropriate local opera- 
tions Rj's. By following the lines depicted in Sec. [Til we get 



the SASA function 

SASA(V, d) = 2Erf 3 



V2d 


^1 + Erf 


y/2d 






^ _ 



(20) 



where perfectly efficient homodyne detectors have been as- 
sumed, for the sake of clarity. Our analysis can be straight- 
forwardly extended to the case of inefficient detectors simply 
by performing the replacements for d and V discussed for the 
case of a tripartite GHZ -like ETS state. Fig.|7]shows the trend 
of Eq. (|20| > against V and d for a considerable range of val- 
ues. For V ~ 10 3 , only moderate values of d are sufficient to 
violate the local realistic bound, although maximum violation 
occurs only at d > 50. Hence our linear cluster ETS maxi- 




SASA(V,d) 



FIG. 7: Violation of SASA inequality by a quadripartite cluster-like 
ETS under effective rotations and homodyne detection. The SASA 
function is plotted against V and d. The local realistic bound is 
shown by the plane at SASA = 2. 
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mally violates this inequality for any value of V by choosing 
d sufficiently large. 

The fact that one of the parties in SASA has to use 
only a single measurement setting (the local observable for 
mode 2 has to be 1 in two cases out of four) marks the 
fundamental difference between such an inequality and any 
Mermin/Svetlichny-like argument, where at least two mea- 
surement settings are required for each of the parties involved. 
Indeed, SASA inequality does not fall into the class of gener- 
alised Mermin inequalities, such as those proposed by Collins 
et al. [20], even more given that it stems from the use of 
stabilizer operators for the cluster state at hand. It is there- 
fore interesting to test which is the behavior of our four- 
mode cluster-like ETS against Mermin-like inequalities. For 
this task, we consider the works by Werner- Wolf 113211 and 
Zukowski-Brukner 113 311 , who independently proposed equiv- 
alent and general inequalities for an arbitrary number of parti- 
cles, using two measurement settings per party. Here we con- 
centrate on the formulation by Zukowski and Brukner, whose 
inequality for spin- 1/2 particles takes the form 

\WWZB\ = \C(abcd)+C(abcd 1 )+C(abc 1 d) - C(abcidi) 

+ C(abicd)—C(abicdi)—C(abiCid)—C(abiCidi) 

+ C(aibcd)—C(aibcdi)—C(aibcid)—C(aibcidi) 

—C(aibicd)—C(aibicdi)—C(aibiCid)+C(aibiCidi)\ < 4 

(21) 

with WWZB indicating the Werner- Wolf-Zukowski-Brukner 
parameter. By using the spin- 1/2 linear cluster state \Ci/ 2 ), 
it is seen that the local observable for particle j needed to build 
up the correlations all have the form 

*to>=V5 -i)> (22) 

so that k = (<pi,4>x) and analogous definitions for the re- 
maining set of parameters. Numerical maximization of the 
left-hand side of Eq. (fJTJ over such sets leads to the viola- 
tion of WWZB inequality by a factor y/2 for the set of angles 
<f>j = 37r/16, (f)'j — ll7r/16, Vj. We retain such values for our 
test of non-locality involving the cluster-like ETS and calcu- 
late the appropriate correlators by means of our dichotomized 
homodyne measurements. The corresponding WWZB func- 
tion is too lengthy to be given here and we thus rely on Fig. [8] 
for a comprehensive account of our results, which show that 
violation up to the maximum allowed to the class of states 
addressed here is achievable under the same general working 
conditions valid in the case of the SASA inequality. 



C. Alternative generation scheme for linear cluster-like ETS 

The protocol for state generation suggested in the previous 
Subsection is not the only possibility one has to generate lin- 
ear cluster-like ETS. Indeed, along the lines of the scheme for 
tripartite GHZ-like ETS proposed in Fig. [U(b), one can de- 
sign situations where a suitably off-line prepared non-classical 
state of two modes is made available for interference at a set 








FIG. 8: Violation of WWZB inequality by a four party linear cluster- 
like ETS under effective rotations and homodyne detection. The lo- 
cal realistic bound is 4 and is indicated by the horizontal plane. For 
any V, a correspondingly large value of d guarantees a violation of 
the WWZB inequality by a factor y/2, in agreement with the spin- 
1/2 case. 



of beam splitters. This is the case for the procedure depicted 
in Fig. [6] (b), where two thermal displaced states are fed into 
a cross-Kerr medium of large third-order non-linear rate x- A 
similar scheme has previously been put forward by Ba An and 
Hoa i34ll . The inter-mode coupling model within the medium 
thus reads H x (z) = hxd\d\d\d^,. By assuming an interaction 
time equal to 7r/x and using the results of Ref. l35ll one finds 
that the output modes are in state 

Pi> = J d 2 ad 2 fiP^(V,d)Pf(V,d) |V) 13 (V| (23) 

with 1^)13 = \ a )\ \ s ^)z + | — a )\ \ s ~p)z- Interference of such 
a state with two vacuum-state modes at 50 : 50 beam splitters, 



SASA(V,d) 




FIG. 9: Comparison between the SASA function associated with 
states p , . (dashed line) and a, . (solid line) for V = 1, 5 

' cluster v ' 1 cluster y ' 

and 10. 
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together with local phase-shifting operations (not shown in 
Fig. 0(b)), leave us with the state 

P { cLter= J d 2 ad 2 f3P*(V,d)P$(V,d) Pc (a',p') (24) 

with Pc (a,/3) = \C(a',a',j3',j3')) 123 4(C(a',a',l3',l3')\ and 
a' = a /V2, P' = /?/ a/2- The crucial difference between 

Pollster and Pollster lies on the lack of independence between 
mode 1 and 3, as well as 2 and 4, which share the same am- 
plitudes. This is the reason behind some difficulties in the full 
analytical treatment of the problem at hand, similar to those 
experienced in treating the situation described in Fig. [TJ (b). 
We thus rely on a numerical analysis. Given that SASA in- 
equality is optimal for linear cluster states, we concentrate on 
this case and compare the performance of the SASA function 
associated with p^r uster and the behavior of Eq. ( f20l > revealed 
in Fig. [7] Features qualitatively similar to those shown in the 
tripartite case are evident in Fig. |9j where the beam splitter- 
based approach is shown to be more advantageous as it allows 
for the violation of SASA inequality at slightly smaller values 
ofd. 



IV. CONCLUSIONS 

We have assessed the problem of revealing genuine multi- 
partite non-locality in many-mode entangled states based on 
chaotic local components. We have identified in Svetlichny 
inequality and stabilizer-based arguments the crucial tools for 
our study, which has been performed by considering relevant 
instances of multipartite entangled states. In particular, we 
have considered tripartite GHZ and W-like ETS as well as 
quadripartite linear cluster states, revealing that the non-local 
properties of the class of states at hand can be faithfully re- 
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